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ABSTRACT 

The X-ray source CXOXBJ142607.6+353351 (CXOJ1426+35), which was identified in a 172 ks 
Chandra image in the Bootes field, shows double-peaked rest-frame optical/UV emission lines, sepa- 
rated by 0.69" (5.5 kpc) in the spatial dimension and by 690 km in the velocity dimension. The 
high excitation lines and emission line ratios indicate both systems are ionized by an AGN contin- 
uum, and the double-peaked profile resembles that of candidate dual AGN. At a redshift oi z = 1.175, 
this source is the highest redshift candidate dual AGN yet identified. However, many sources have 
similar emission line profiles for which other interpretations are favored. We have analyzed the sub- 
stantial archival data available in this field, as well as acquired near-infrared (NIR) adaptive optics 
(AO) imaging and NIR slit spectroscopy. The X-ray spectrum is hard, implying a column density 
of several 10^^ cm~^. Though heavily obscured, the source is also one of the brightest in the field, 
with an absorption-corrected 2-10 keV luminosity of ~ IC*^ erg s~^. Outflows driven by an accretion 
disk may produce the double-peaked lines if the central engine accretes near the Eddington limit. 
However, we may be seeing the narrow line regions of two AGN following a galactic merger. While 
the AO image reveals only a single source, a second AGN would easily be obscured by the significant 
extinction inferred from the X-ray data. Understanding the physical processes producing the com- 
plex emission line profiles seen in CXOJ1426-I-35 and related sources is important for interpreting the 
growing population of dual AGN candidates. 

Subject headings: galaxies: active - galaxies: nuclei - quasars: emission lines - galaxies: individual 
(CXOXBJ142607.6-I-353351) 
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1. INTRODUCTION 

It is now generally believed that galaxies evolve in a hi- 
erarchical fashion, with galaxies growing through merg- 
ers. Furthermore, evidence suggests that most, if not 
all, galaxies con tain a supermassive b lack hole (SMBH) 
at their center (jRichstone et al.l 119981) . Therefore, in a 
galactic merger the force of dynamical friction will act 
to bring each SMBH toward the center of the new mass 
distribution, first forming a dual SMBH (separations of 
several kpc), then a gravitationally bound system (sub- 
pc separations) and eventually merging to form a sin- 
gle, more massive SMBH. The amount of time spent in 
each stage is dependent upon the efiiciency of the pro- 
cesses carrying away angular momentum from the system 
(|Begelman et al.lll98gi . At large separations (> 10 — 100 
pc), the primary force acting on the SMBHs is dynam- 
ical friction from the surrounding gas and dust. As 
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the SMBHs form a binary system (i.e. harden), ejec- 
tion of stars becomes most important. Finally at close 
separations (< 10~^ pc), gravitational radiation will 
become the most important m echanism of energy loss 
(jMilosavlievic fc Merr'itt ''20031. 

Numerical simulations show that in a wet merger (i.e., 
where the merging galaxies possess a significant fraction 
of gas and dust), tidal interactions are effective at fun- 
neling the gas and dust toward the center of the merg- 
ing system where the two SMBHs are located, provid- 
ing fuel for enhanced accretion (Hcrnquist 1989). This 
process is thought to be a possible mechanism for initi- 
ating an active galactic nu cleus (AGN) or quasar phase 
(jHopkins et al.l [20051 12008D . and the existence of a dual 
AGN (two actively accreting SMBHs within the same 
galaxy) is a natural prediction of this merger fueling 
model. Therefore, identifying pairs of SMBHs in the pro- 
cess of merging is extremely important for probing the 
link between AGN activity and galaxy evolution. 

Unfortunately, observational evidence of such systems 
has proven very difficult to obtain. Large samples 
of binary quasars exist as evidence for quasar c l uster- 
ing and mu t ual t r iggeri n g (iHennawi et al.l 120061 120101 : 
IMvers et al.l [20071 [200l [Foreman et al.l 120091) . How- 
ever, these are typically at large projected separa- 
tions (> 10 kpc) and thus generally only represent 
the earliest stage of binary SMBH evolution. Most 
confirmed cases of dual SMBHs at smaller separa- 
tions (between 10 kpc and 1 kpc) have been identi- 
fied in local ultraluminous infrared gala xies (ULIRGs) 
as pairs of X-ray point sources (Komos sa et al.l 120031: 
iGuainazzi et al.l[2005t [Hudson et al.M2006l : iBianchi et al.l 



2 



Barrows et al. 



[20081: iPiconcelli et al.llMoHKoss et al.ll2011[ ). with much 
of the optical eniission obscured by the heavy internal ex- 
tinction common in ULIRGs. In some cases the obscu- 
ration is so heavy that only one nucleus is visible, thus 
highlighting the difficulty that dusty, merging systems 
present to identification of dual AGN. Additionally, the 
close (7 pc separation) binary SMBH system in the radio 
galaxy 0402-1-379 was identified fro m Very L ong Base- 
line Array (VLB A) observations (Ro driguez et al. 2006). 
While the excellent spatial resolution achievable in the 
radio regime makes the VLBA capable of detecting close 
binaries, this technique is time consuming and not suited 
to identifying a large sample of sources. 

In the past few years, several interesting individ- 
ual examples of dual/binary SMBH candidates have 
been identified from single-epoch optical spectra of 
AGN or quasars showing double-peaked emission lines. 
The Sloan Digital Sky Survey (SDSS) spectral archive 
(|Abazaiian et al.|l2009[ ) has been a fruitful database for 
such identifications. However, these objects do not rep- 
resent a homogeneous class, and thus the physical inter- 
pretation of many of them, based on single-epoch fiber 
spectra alone, remains ambiguous. 

For example, the sources SPSS J15363 6.22-I-044127.0 
(J1536+0441, IBoroson fc Laue^ 120091) and SDSS 
J10504 1.35-H345631.3 (J1050-F3456, IShields et al.1 
l2009b() show double-peaked, broad Balmer emission 
lines (Ha, H/3 and H7), which resemble two broad 
line regions (BLRs) with significant velocity differences 
(> 2000 km s~^). However, these objects might in- 
stead be members of the sub-class of AGN/quasars 
called 'double-peaked emitters ' (also known as 'disk- 
emitters') where the double peaks of the Balmer lines 
are attributed to emission from a single accretion disk. 
Line emission contribution from an accretion disk is 
suggested by the extended wings in the Balmer emission 
lines (particularly evident in Ha, e.g. ' Chornock et al.l 
[20T0I) which is a characteristic disk profile feature. On 
the other hand, these sources are extreme examples 
of 'disk emitters,' and it has been shown that such 
strongly double-peaked Balmer lines in quasars may 
represent cases in which the presence of a secondary 
SMBH is effective at enhanc ing the double-peaked 
profile produced by the disk (jTang fc Grindlavl 120091 : 
] Mm. One su ch object, 4C-F22.25 
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identified bv iDecarli et al.l ()2010l ) as a candidate binary 
SMBH, shows very high (8700 km s~^) velocity offsets 
between the two peaks of the Balmer emission lines; 
it is therefore unlikely to be a 'disk-emitter' as the 
velocity difference is several thousand km s~^ larger 
th an that of any known similar source (see the sample 
of iStrateva et al.l l2003| ). These results demonstrate 
the diversity of double-peaked broad line profiles and 
emphasize that in a close binary SMBH there are likely 
to be strong interactions between the SMBHs and the 
accretion disks. 

A growing number of objects have been identified with 
narrow forbidden emission lines that are either offset 
from the systemic velocity or double-peaked in single- 
epoch spectra, many of which have been found through 
systematic searches for candidate du al AGN, both in 
the D EEP2 Galax y Redshift Survey (iComerford et al.l 
2009al) and SPSS (IZhou et aLI [2003 IWang et al.l 120091: 



jects have been identified as dual AGN candidates be- 
cause they show s t rongly double-peaked [OIII]5007 line s 
(|Zhou et all I200I IGerke et al.l[2007t IPeng et aLllMll) . 
or double-peaks in all of the narrow lines, both per- 
mitted and fo rbidden, e.g. SPSS 9 2712.65+294344.0 
(J0927 +2943. iKomossa et al.l l2008bl : iBogdanovic et al.l 
I2009bl: iPotti et al.ll2009[ ). These double-peaks indicate 
two sets of narrow emission line regions (NLRs), and 
are suggestive of dual AGN with separations up to sev- 
eral kpc. These systems are excellent sources for study 
because, at the implied kpc-scale separations, the two 
putative AGN can be resolved in two-dimensional spec- 
troscopy and imaging. 

However, double-peaked narrow lines have been seen 
in a number of AGN and quasars where other phys- 
ical explanations have been favored. For example, 
NLRs are thought to have complex gas kinematics which 
may include radially flowing material capable of pro- 
ducing offset or doub le-peaked narrow emission lines 
(|Heckman et al.l I1981D . In particular, powerful radio 
galaxies often show complex, spatially extended NLRs 
with multiple components aligned along the radio axis 
(see reviews of rad io galaxies in IMcCarthvl 119931 and 
iMilev fc Pe BreucM 2008). Models for t his alignment in- 
clude reflection of the AGN emiss ion (iTa dhunter et all 
1198 ^. jet-induced star formation (iPe Yo ung 1981]), in- 
verse Compton scattering of cosmic microwave back- 
ground photons of f the relativistic electrons in the ra- 
dio jet (jPalvl [199^ ■ material out- flowi ng from an acce- 
tion disk or entrained in a radio jet (|Holt et al.l I2003L 
l2008t IKomossa et al.ll2b08a|) . as well as a combination 
o f shocks and photo ionization by the central engine 
(iiVIaxfield et al.l 120021) . There is an extensive litera- 
ture on the spatially extended, k inematically complex 
NLRs of powerful radio galaxi e s (Ivan Oiik et al.l ^99^ 
VUlar-Martin et al. 2000l 120071 I2OIOI: iZirm ~ al. 2005; 
Nesvadba et al.l 12006 ) . Though radio galaxy NLRs of- 
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ten appear vaguely reminiscent of dual AGN candidates, 
multiple AGN are rarely suggested as possible interpreta- 
tions of the complex radio galaxy NLRs. Finally, simple 
geometrical interpretations have been invoked to explain 
double-peaked narrow lines in AGN, such as a rotating, 
disk-shaped NLR or a single AGN illuminating two NLRs 
(|Xu fc Komossall20iDa . 

Peveloping a clear understanding of the physical pro- 
cess producing these often complex emission line pro- 
files is important to both the search for dual AGN and 
to further understand quasar/ AGN emission line regions 
in general. Since it is clear that a number of physical 
processes other than a pair of AGN can explain double- 
peaked emission lines (broad and/or narrow), more infor- 
mation is necessary to entirely reject or confirm any one 
model. From single-epoch spectra, emission line ratios 
may be used as diagnostic tools to constrain the nature 
of the ionizing continuum producing each set of emis- 
sion lines. Other observational tools include spatially re- 
solved spectroscopy, synoptic spectroscopic observations, 
and high resolution imaging. 

In this paper we discuss the X-ray source GX- 
OXBJ142607.6-f 353351 (from here on CXOJ1426-h35) 
which was originally identified as an X-ray source in 
the Bootes field (Wang et al. 20Q1). In follow-up slit- 
spectroscopy, CXOJ14264-35 showed double rest-frame 
optical/UV emission lines, separated spatially by 0.69", 
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and in velocity-space by 690 km s ^, as measured from 
the [OIII]5007 line. We discuss CXOJ1426+35 as an ex- 
ample of a candidate dual AGN, and place its spectral 
and morphological properties within the context of other 
candidate dual/binary SMBHs. We examine a variety 
of possible explanations for the double-peaked emission 
lines. Interestingly, CXOJ1426-f35 is at z = 1.175, the 
highest redshift yet for a candidate dual AGN, where 
previously candidates were identified out to z ss 0.8. At 
this redshift, we have access to rest-frame UV emission 
lines which is unique among the current sample of can- 
didate dual AGN in the literature. Section [2] discusses 
the known properties of CXOJ1426-I-35 and Section [3] 
discusses various physical interpretations of the data. 
Section [4] takes CXOJ1426+35 as a pedagogical exam- 
ple to describe a number of methods by which candidate 
dual/binary SMBHs may be confirmed or rejected. Sec- 
tion [5] contains the conclusions. Throughout the paper 
we assume cosmological parameters of Hq — 70 km 
kpc-\ flM = 0.30 and Qa = 0.70. 

2. PROPERTIES OF THE SOURCE 

In an attempt to understand the nature of the com- 
plex line emission of CXOJ1426-I-35, we first summarize 
in this section the known properties of the source. The 
multi-wavelength data available for CXOJ1426-I-35 in- 
clude both archival data and newly obtained data from 
the Keck II telescope in Hawaii. 

2.1. Archival Data 

The photometric properties of CXOJ1426-f35 from 
archival observations are summarized in Table [T] The 
mid-IR through near-ultraviolet (NUV) colors are dis- 
played as a broad-band spectral energy distribution with 
best-fitting spectral templates in Section [231 We begin 
by describing the archival observations. 

2.1.1. Chandra X-ray Data 

CXOJ1426-I-35 was detected with the A dvanced CCD 
Imaging Spectrometer fACIS-I: IGarmire" et al. 2003) on 
board the Chandra X-Ray Observatory (CXO) in a 172 
ks im age of the Large A rea Lyman Alpha (LALA) Bootes 
field (|Wang et al.ll2004l) . The total image was composed 
of two expos ures take n on A pril 16-17 and June 9 of 
2002. In Wa ng et all ()2004f) . the reported number of 
background-subtracted soft (0.5 - 2.0 keV) and hard (2.0- 
7.0 keV) counts, determined by 'wavdetect,' are 2A.3'!^f l 
and 240.8l^]' o, respectively, which yields a hardness ra- 
tio (defined as HR = {H - S)/iH + S), where H and S 
are the number of the hard and soft counts, respectively) 
of 0.82. Assuming Galactic absorption and a power law 
spectrum, the hardness ratio corresponds to a photon in- 
dex of r = —1.3, which is significantly harder than the 
typical observed F = 1.4 of AGN when not corrected for 
absorption. Interestingly, out of 168 sources detected in 
the field by I Wang et al.l (|2004f ). only thirteen have larger 
0.5-10 keV fluxes and only seven are more obscured (i.e., 
have higher hardness ratios); no brighter sources have 
higher hardness ratios. Thus, CXOJ1426-I-35 is unique 
in that it is bot h relatively bright and heavily obscured 
(i.e.. Figure 2 of lWang et 311120041) . 

Unfortunately, the off-axis detection of CXOJ1426+35 
means that it is unresolved by CXO. However, the high 



number of hard counts allowed us to obtain a reason- 
able spectrum. To do so, we re-reduced the data us- 
ing the CIAO software, following the standard proce- 
dure for extracting source and background regions, and 
creating response files for spectral analysis. The two ob- 
servations were processed independently. In each event 
file, we extracted a source region 30 pixels (15") in di- 
ameter centered on the centroid position of the source 
and a background region 50 pixels (25") in diameter ad- 
jacent to the source region. The source regions, back- 
ground regions and response files were used in XSPEC 
to perform spectral modeling. For each event file, the 
spectra were re-binned so that each bin had at least 15 
counts. When fitting the models, the same model was fit 
to the spectrum from each observation simultaneously 
to place tighter constraints on the parameters. All bad 
channels were ignored, and the first channel of each data 
set (where there are few counts) was also ignored. A sim- 
ple absorbed power law was fit with the redshift frozen 
at z = 1.175 (see Section [2.2.11 for the redshift deter- 
mination). The spectrum, with the model overlaid and 
residuals underneath, is shown in Figure [T] The resid- 
uals around ~ 3 keV are suggestive of rest-frame 6.4 
keV Fe Ka emission, though there are too few counts 
to obtain a good line fit. The best fit parameters are 
an intrinsic power-law slope F = 1.961^ 57 ^^'^ column 
density nn = 5.90l^'8^ x 10^^ cm^^. This column den- 
sity implies a highly obscured source. From the best fit 
model the observed 2-10 keV flux is 5.44l°;^^ x 10"" erg 

s-i cm-2 or 2.8lt;5ji^ x lO'^ ^Jy (Table |^ which cor- 
responds to an absorption-corrected intrinsic luminosity 
of L2-iokoV = 9.7l:2;3 X 10' ^^ erg s~^ From Lz -iokoV, 
the bolometric correction of iHopkins et al.l (|2007l ) yields 
Lbol « 9 X 10*6 erg s-\ 

2.1.2. UV Imaging 

CXO J 1426-^35 was detected by the NU V detector on 
the G alaxy Evolution Explorer {GALEX; IMartin et al.l 
I2003D which has an observed bandwidth of 1771-2831 
A and effective wavelength of 2271 A. The 5" FWHM 
resolution of the NUV detector on GALEX is unable to 
resolve CXOJ1426-I-35. Therefore, for photometry we 
have used the flux measurement from the 3.8" radius 
aperture (to most closely match the 4" radius aperture 
for the IRAC data; Section I2.1.4p . We then corrected 
to a total magnitude using the correction factor from 
iMorrissev et al.l (|2007D and converted to a Vega magni- 
tude calculat e d bas ed on the model atmosphere of Vega 
from ' Kurucd (|1993l ). The NUV photometry is listed in 
Table [TJ CXO J 1426-^35 was not detected by the far- 
ultraviolet (FUV) detector on GALEX. 

2.1.3. Optical/Near-IR Imaging 

CXOJ1426+35 is well-detected by the SDSS (Figure 
[2]). In Table [1] we list the model magnitudes, calcu- 
lated by the SDSS pipehne using the best fitting func- 
tion (either a deVaucouleurs or an exponential profile 
- in this case an exponential profile) in the r band 
and applying that model to the other bands. We con- 
ver ted the SDSS magn itudes to AB magnitudes foUow- 
ing lBohhn etal\ (pOOll ): ABoSset « 0.04, 0.0, 0.0, 0.0, and 
—0.02 magnitudes for u, r, i, and z, respectively, where 
-^^AB = AfsDss^ABoftsot- We then converted to Vega 
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magnitudes ba sed on the model atmosphere of Vega from 
iKurucd ()1993[ ). Figure [2] shows the stacked g + r + i 
SDSS image of CXOJ1426+35, reveaUng a clearly ex- 
tended host galaxy. 

Note that the 8.5 deg^ Bootes field has been deeply im- 
aged at optical wavelen gths by the NOAO D eep, Wide- 
Field Survey fNDWFS: Uannuzi et al.|[2000l) . However, 
the bright star 25" East of our target (see Figure [2]), 
while useful for the AO observations discussed in Section 
12.2.21 severely contaminates the NDWFS photometry for 
CXOJ1426-I-35 due to a bright diffraction spike. 

The Bootes field has also been deeply imaged at NIR 
wavelengths (J, H, and Kg) using the NOA O Extremely 
Wide-Field Infrared Imager (NEWFIRM; Prob st et all 
12004 ,2008) on the 4-m Mayall Telescope at Kitt Peak 
National Observatory. Those data were obtained as 
part of the Infrared Bootes Imaging Survey, an NOAO 
Large Survey Program (Gonzalez et al., in prep.). 
CXOJ1426-h35 is detected in all bands. Tableffllists the 
corresponding photometry in both Vega and AB magni- 
tudes, where the latt er assume the Vega-AB offsets from 
iBlanton et all (pOOSi ). 

2.1.4. Mid-IR Imaging 

CXO J 1426-1-35 has mid-IR data availa ble from the 
SpUzer Deep, Wide-Field Survey (SDWFS; lAshbv et al.l 
l2009f l. a four-epoch infrared survey of 10 deg^ in the 
Bootes field usin g the Infrared Array Camera (IRAC; 
iFazio et aLll20Q4[ ) on the Spitzer Space Telescope. Mag- 
nitudes are provided for channels 1, 2, 3 and 4, which 
correspond to bandpasses with effective wavelengths of 
3.6, 4.5, 5.8 and 8.0 fim, respectively. The photome- 
try listed in Table ffl assumes the position measured in 
the 4.5 /xm image, and uses 4.0" radius aperture pho- 
tometry, corrected to total magnitudes, and converted 
from Vega magnitudes to AB magnitudes using constants 
listed in the IRAC Data Handbook. The mid-IR colors 
are AGN-like, falling within the empirical separation of 
active galaxies fro m Galactic stars and normal galaxies 
(jStern et al.ll2005[ ). though it does lie on the red end of 
the overall AGN distribution. 

2.1.5. Radio Imaging 

CXOJ1426-f35 is undetected in the deep, 1.4 GHz ob- 
servation of the Bootes field by the W esterbork Synthesis 
Radio Telescope ()de Vries et al.|[200^ . These data reach 
a median noise level of 140/iJy {5a), and at this depth, no 
sources are detected within 14.6' of CXOJ1426-h35. For 
a source at a redshift of z = 1.175, this detection limit 
corresponds to a radio power of Pi. 4 ghz = 1 x lO^'^ W 
Hz""'^. For comparison, we used th e 1.4 GHz radio p owers 
of the 14 compact radio sources in IHolt et al.l (l2008h (cal- 
culated using the published power law indices) as those 
sources also exhibit narrow emission lines with velocity 
offsets, presumably as a result of interactions between 
the radio jets and the NLR clouds. The range of radio 
powers is P1.4 ghz = 3.8 x 10^^-1.5 x 10^^ W Hz^i, with 
a median of Pi. 4 ghz = 3.8 x 10^^ W Hz~^, meaning this 
sample of compact radio sources would have all been de- 
tected in this observation. Therefore, CXOJ1426-I-35 is 
unlikely to be a similarly strong, compact radio source. 

2.2. New Data and Analysis 

While targeting a mid-IR transient w hich proved to be 
a self-obscured supernova at z 0.2 (jKozlowski et al.l 



[2OT0h . we obtained a deep spectrum of CXOJ1426-I-35 at 
Keck Observatory in March 2010. As described below, 
this intriguing spectrum triggered several additional 
observations. 



2.2.1. Keck/LRIS: Optical Spectroscopy 

Optical spectra of CXOJ1426-I-35 were obtained 
with the L ow Resolution Imaging Spectrometer (LRIS; 
lOke et all [19951 on UT 2010 March 12 using the D680 
dichroic and the 400/3400 grism on the blue side and 
the 400/8500 grating on the red side. The 1.5" slit was 
oriented at P.A.= —21°. The data consisted of three ex- 
posures of 1200 seconds each. For a 1" slit, this instru- 
ment configuration provides a resolution of 6.5 — 7.1 A 
(FWHM) on the blue side of the dichroic and 6.9 A 
(FWHM) on the red side of the dichroic for objects fill- 
ing the slit. We processed the data using standard tech- 
niques within IRAF. The 2D-spectra (Figure[3]) show two 
emission line systems separated both spatially and along 
the dispersion axis. At a P.A. of —21°, the spatial sepa- 
ration is 0.48", as measured from the [NeV]3426 line. To 
model the emission line velocity-offsets and widths, we 
extracted an aperture of 3" in diameter to recover the 
flux from the entire source (both emission line systems 
plus the continuum), and separate apertures centered on 
each component to isolate the flux from the 'blue' and 
'red' systems. Figure S] shows the 3" extraction with the 
emission line identifications. These extractions mimic 
what a SDSS fiber spectrum of the source would look 
like, though the spatially resolved spectrum, shown in 
Figure [21 provides considerable additional information. 

The extracted spectra were corrected for Galactic ex- 
tinction using a color ex cess of E(B — V) = .015 deter- 
mined from themodel of iSchlegel et~all (|1998f ). We used 
SPECFIT (jKrissI 119941 ) to model each emission line as 
two Gaussian profiles (for the 'blue' and 'red' emission 
line systems, respectively) on top of a local continuum. 
Unfortunately, the individual component apertures have 
significant contamination from each other, making it dif- 
ficult to deconstruct the emission line fiux attributable 
to each emission line system. Therefore, those apertures 
were used only to obtain estimates for the observed cen- 
tral wavelengths of each emission line, and the lines were 
fit again in the 3" aperture to determine the spectral 
model parameters and associated errors. For the highly 
blended emission lines, some of the Gaussian parameters 
were fixed to those of less blended lines with similar ion- 
ization potentials. Doublets with rest wavelengths close 
together were modeled as single Gaussians. We also de- 
tected an absorption feature which is consistent with an 
H9 absorption line at a redshift intermediate to the two 
emission line systems. For each free parameter that was 
fit, the quoted errors correspond to la uncertainties, and 
these errors were propagated into any further calcula- 
tions. The line fluxes, FWHMs, redshifts and velocities 
separations are presented in Table 121 

For both the 'blue' and 'red' systems each emission 
line is well-fit by a single Gaussian. As measured from 
[NeV]3426, the redshifts for the 'blue' and 'red' systems 
are z = 1.1715 and z = 1.1773, respectively. The range 
of spectral coverage is 3100 — 5200 A on the blue side 
of the dichroic and 6700 — 9600 A on the red side of 
the dichroic, which corresponds to rest-frame wavelength 
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ranges of 1430 - 2380 A and 3080 - 4410 A and allows 
us access to a range of rest-frame optical and UV emis- 
sion lines that are useful in diagnosing the structure and 
kinematics of the NLR gas from their velocities and line 
widths. If the absorption line is interpreted as H9 absorp- 
tion within the host galaxy, this yields a systemic redshift 
of z = 1.1751. At this redshift, the emission line velocity- 
splitting is 790 km as measured from [NeV]3426. 
With no broad components detected in any of the op- 
tical/UV emission lines, the spectrum of CXOJ1426+35 
resembles that of a Type 2 AGN, which in the unified 
theory of AGN implies that the BLR and central con- 
tinuum source are obscured. This is consistent with the 
heavy obscuration implied by the X-ray data. 

2.2.2. Keck/NIRC2: Adaptive Optics Imaging 

Adaptive optics (AO) imaging of CXOJ1426+35 was 
obtained by the na rrow camera of the Near Infrared 
Camera 2 (NIRC2; Ivan Dam et"all I2004D at the Keck 
II telescope on UT 2010 June 30. These data con- 
sisted of three 180 second data exposures at P.A.= 0.0° 
through the K filter with an effective wavelength of 
Aeff = 2.12 /im which corresponds to a rest- frame wave- 
length of Aoff ~ 1/im. An average dark field image was 
created using 'imcombine' and was subtracted from each 
data image. Each image was flat-fielded and then differ- 
enced from the other science frames. The images were 
shifted using the position of the target and then stacked 
using 'imcombine'. The final stacked image, 10" on a 
side, is shown in Figure [5] and compared to the g -\- r + i 
SDSS image and the NIRSPEC image (see SectionHXS]). 
The image reveals a single, diffuse source. The lack of 
two nuclei in the K band implies that either there is 
no second AGN or alternatively that there is significant 
obscuration around l/im. Significant obscuration would 
also be consistent with the high column density implied 
by modeling of the X-ray data (Section I2.1.ip and with 
the extinction implied by our modeling of the spectral 
energy distribution (Section 12. 3|) . 

2.2.3. Keck/NIRSPEC: Near-IR Imaging and Spectroscopy 

NIR imaging and spectroscopy of CXOJ1426+35 was 
obtain ed using the N ear Infrared Spectrometer (NIR- 
SPEC; IMcLean et "all [1998,) at the Keck II telescope on 
UT 2010 July 18. 

Imaging: The imaging consisted of 11 dithered images 
of 10 seconds each. The filter used was the NIRSPEC-1 
filter which has a wavelength range of 0.95 — 1.12 /im 
and is approximately a Y-band filter. The images were 
flat-fielded and then differenced. The background was 
subtracted using 'background' for each quadrant of each 
image independently and cosmic rays were removed using 
'xzap' in IRAF. The images were stacked with 'imcom- 
bine' using for offsets the brightest star, when available, 
and the target otherwise. At a redshift of z = 1.175, the 
fllter provides rest-frame coverage of 0.44 — 0.52 fim. The 
flnal, combined image is presented in Figure [S] and reveals 
that the host galaxy dominates at optical wavelengths 
and consists of an elongated and "lumpy" structure that 
hints at a disturbed morphology. 

Spectroscopy: We obtained NIRSPEC spectroscopy in 
low- resolution mode using a 0.7" wide slit placed at a po- 
sition angle of P.A.= —54.1° so that the spatial axis was 
parallel to the major axis measured in the SDSS image. 



For a 0.38" slit, the spectrograph provides a resolving 
power of ~ 2, 200. At two positions, A and B, along 
the slit, four 300 second exposures were taken in the or- 
der ABBA. The same fiher as for imaging, NIRSPEC- 
1, was used. The spectra were differenced, cleaned of 
cosmic rays using the IRAF task 'szap' and background 
subtracted using 'background'. The spectra were shifted 
and stacked using 'imcombine'. 

The range of spectral coverage allows us access to the 
H/3, [OIII]4959 and [OIII]5007 lines, which clearly show 
two emission line components (Figure [S]). As measured 
from [OIII]5007, the separation between the two compo- 
nents along the spatial axis (parallel to the major axis of 
the galaxy in the NIRSPEC image) is 0.69". Aperture 
extractions and spectral modeling proceeded as with the 
LRIS spectrum. Three apertures were extracted that in- 
clude a region of 3" diameter encompassing the entire 
source, and separate apertures for the each of the 'blue' 
and 'red' systems. The extracted spectrum was corrected 
for Galactic extinction, and spectral modeling was per- 
formed on the 3" aperture with 'SPECFIT' with initial 
guesses for the line central wavelengths determined from 
the individual apertures for each component. The 'red' 
II/3 component is very weak and therefore we have fixed 
its FWHM to that of the 'red' [OIII]5007 component and 
adopted the measured fiux as an upper limit. Line fiuxes, 
FWHMs, redshifts and emission line velocity-splittings 
are listed in Table [2j 

The emission lines for both the 'blue' and 'red' sys- 
tems are each fit by a single Gaussian (Figure [7]) • The 
[OIII]5007 line has the highest S/N and the wavelength 
solutions for the Gaussian peaks provide a line of sight 
velocity separation of 690 km s^^. The ratios between 
il/3 and H7 are consistent within their errors with case B 
recombination theory for a reasonable NLR temp erature 
range of 2,500-20,000 K (jOsterbrock fc Ferlandl [20061. 
Considering the significant error of the II7 flux measure- 
ment, we do not find any measurable Balmer decrement 
in the narrow emission line systems. The ratio between 
[OIII]5007 and H/3 is highly useful in diagnosing the na- 
ture of the ioni zing continuum throu gh its placement on a 
BPT diagram (iBaldwin et al.lll981f ). The line fiuxes pro- 
vide [OIII]5007/H/3 ratios of 8.9±3.2 for the 'blue' system 
and a lower limit of 11.5 for the 'red' system. Though 
we do not have access to Ha and can not uniquely place 
the en iission line system s on a typical diag nostic dia- 
gram (iKauffmann et al.l 120031 : iKewlev et a"lTl2006[ ). the 
[OIII]5007/H/3 r atios put both systeni s above the demar- 
cation defined in IKewlev et al.l ()2001t ). clearly indicative 
of AGN and consistent with the high ionization lines ob- 
served in the optical spectrum. 

2.3. Modeling of the Spectral Energy Distribution 

Motivated by the high level of AGN obscuration im- 
plied from the X-ray HR and by the SDSS and NIR- 
SPEC images, which suggest that the galaxy dominates 
at optical wavelengths, we have attempted to separate 
the SED into galaxy and AGN components to deter- 
mine their relative contributions. This was done to de- 
termine if the AGN extinction i nferred from the S ED is 
consistent with the X-ray HR ^Hickox et al.l [20071) and, 
furthermore, if the SED is potentially consistent with 
two AGN. We modeled the broad-band photometry of 
CXOJ142 6-K35 using the SE D templates for AGNs and 
galaxies of lAssef et al.l ()2010f ) . The four empirically con- 
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structed SED templates that form this basis extend from 
0.03 to 30 /im and include an old stellar population (re- 
ferred to as E), an intermediate star forming popula- 
tion (Sbc), a strongly star- forming populatio n (Im), and 
an AG N component. We refer the reader to lAssef et al.l 
(|2010H for d etails on the SEP templates. The standard 
approach of lAssef et al.l ()2010l ) is to fit the broad-band 
photometry with a non-negative linear combination of 
these four templates, allowing for additional dust extinc- 
tion for the AGN component (and only the AGN compo- 
nent). However, this approach does not allow for a satis- 
factory fit to the broad-band photometry (x^ = 185; Fig. 
[S]a). Led by the possibility of CXOJ1426+35 being a dual 
AGN, we have modified our analysis and attempted to 
model the SED with two independent AGN components, 
each with a different amount of extinction, in addition to 
the host components. This approach yields a much better 
fit (x^ — 54), as shown in Figure[8j). An unreddened low 
luminosity AGN component is used to model the UV ex- 
cess over the galaxy-host SED, while another AGN com- 
ponent, intrinsically 10 times more luminous and with 
large amounts of obscuration (E(B ~V) = 3.44), is used 
to model the mid-IR emission. This fit, while much bet- 
ter than that obt ained with the standard approach of 
lAssef et al.l (j2010H , is still not completely satisfactory; 
the «— band, which includes redshifted [OII]3727 emis- 
sion, is underpredicted by ^ 3.5 cr, and the IRAC 8.0/xm 
band is underpredicted by ~ 4 cr. This does, however, 
lend credibility to the hypothesis that CXOJ1426-I-35 is a 
dual AGN. The SED modeling is compatible with the X- 
ray HR - the color excess implies rif f » 2 x 10^^ cm~^, as- 
summg an SMC-like gas/dust ratio (|Maiolino et al.l200H ) 
- and with the single source detected in the AO image, 
as the more luminous but highly reddened AGN compo- 
nent would dominate the X-ray emission but would be 
virtually invisible at rest-fr ame 1 /nm (Aijx m = 4.66 mag 
from the extinction law of (iCardelli et al.llTQSgf ) and for 
Ry =3.1.). The optical broad lines of the less luminous, 
unreddened AGN (with ~ 10% of the host luminosity ) 
would be too faint to be detectable (|Hopkins et al.ll2009l ). 
consistent with our observations. We also find no signif- 
icant evidence for a broad component in the UV lines, 
though in this model a broad component in CIV1549 
would possibly be detectable with sufficient signal qual- 
ity. The relative luminosities of the two AGN compo- 
nents in the dual AGN model are consistent within the 
factor of 4 uncertainty of the bolometric luminosities, 
5 X 10**^ erg s~^ and 1.5 x 10*^ erg s^^ for the 'blue' and 
'red' components, respectively, fr om the -£/[oiiii b olomet- 
ric correction for Type 2 AGN in lLamastra et all (|2009( ). 

We note that an alternative approach would be to asso- 
ciate the UV excess with a very young stellar population 
instead of a second, lower luminosity AGN. Such a sce- 
nario would imply a star-formation rate of approximately 
40 Mq yr~^ us ing the total UV lurn inosity at 2150 A and 
the relation of lMadau eFaLl ()1998D for a Salpeter IMF, 
and a stellar mass of approximatel y 10^° us ing the 
g—r color and M/Lk relation of Be ll et al.l (|2003[ ) for the 
host. A specific star- formation rate of ^ 4 x 10~^ yr~^ 
is close to the high end of the distribution at the redshift 
of CXOJ1426-H35 fsee iCooper eralll2008[) but it is not 
implausibly large, making this modeling of the SED a 
viable alternative to the dual AGN scenario. 



3. INTERPRETATION 

In this section, we discuss six physical scenarios which 
may account for the unusual double emission line peaks 
observed in the restframe optical/UV lines, and the mor- 
phology as evident from the images, particularly the lack 
of a double nucleus in the NIRC2-A0 image. These sce- 
narios are a chance superposition, a rotating accretion 
disk, jet/cloud interactions, an unusual NLR geometry, 
accretion disk outflows, and a dual AGN. Similar expla- 
nations have been invoked for other dual AGN candidates 
and sources with double-peaked or offset emission lines 
(in particular, see th e discussions of Gerke et al. 200"^, 
IKomossa et"al]|2008al and lXu fc Komossa 2009). though 
we elaborate upon them here. We discuss the applicabil- 
ity of each scenario to CXOJ1426-I-35, concluding that a 
dual AGN is a viable explanation. 

3.1. Superposition 

A chance superposition of two unrelated sources has 
occasionally been posited as a possible interpretation of 
the multiple narrow emission lines seen in dual AGN can- 
didates. Ignoring the fact that, in this case, the inter- 
pretation is rather forced since two systems with this 
little separation both physically and in velocity space 
are unlikely to be unrelated, we can use the surface 
density of AGN to estimate the probability of such a 
chance superposition. Using the SDSS quasar catalogue 
(jSchneider et al.ll200 7') to determine the surface density 
of luminous (brighter than Mi = —22.0) quasars at 
1.165 < z < 1.185 (e.g., a redshift bin equivalent to 
nearly one hundred times the velocity difference between 
the two systems), we find a probability of ~ 3 x 10~* 
that CXOJ1426-f35 is due to a chance superposition of 
field galaxies. Though the SDSS quasars are luminous 
Type 1 quasars and do not include obscured sources 
such as CXOJ1426-I-35, this small probability strongly 
argues against a superposition. However, the probabil- 
ity would increase were CXOJ1426+35 a member of a 
cluster (e.g. iDotti & Ruszkowski 2010). Using the clus- 
ter surface den sity at the high end of the surface density 
distribution in iDressler et al.l ()1997f) , we find that the 
probability increases by a factor of ~ 15, implying that 
even in a cluster this scenario is highly improbable. 

3.2. Jet/Cloud Interaction 

Radio jets produced by the central active SMBH 
are capable of introducing a variety of kinematics into 
the NLR. Well-coUimated jets can produce asymmetric 
and/or blueshifted narrow emission line profiles, and in- 
deed jet/cloud interactions in the inner NLR have been 
used to explain th e [QUI] blueshifts seen in some com- 
pact radio sources (IHolt et al.ll200a 12001 iRi fc StocktonI 
|2009() . If the redshifted jet is not obscured, this scenario 
can explain the double-peaked narrow emission lines of 
CXOJ1426+35, and this was shown to be the cas e for a 
previous dual AGN candidate (jRosario "etaDIMoh . Nev- 
ertheless, some questions remain as to the capability of 
radio plasma to entrain the bulk of the NLR so that 
the velocity-shifted lines dominate in the emergent spec- 
trum, rather than introd uce inst abilities that fragment 
the clouds. For example, iXu fc K omossa (2009) discuss 
for another dual AGN candidate the scenario in which 
only a single cloud interacts locally with each side of the 
jet, but find that the [OIII]5007 line luminosity is not 
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consistent with just two clouds, and we arrive at the same 
conclusion for CXOJ1426+35 given its high [OIII]5007 
luminosity. 

However, an even tighter constraint on the possibility 
of a jet/cloud interaction in CXOJ1426+35 is derived 
from the radio power corresponding to the radio detec- 
ti on limit (Section 12. 1.5|) . The sample of radio sources 
in iHolt et ahl (12008 1 would have all been well-detected. 
Though these are only 14 sources and do not represent 
all radio galaxies, they are useful for comparison as they 
show blueshifted [OIII]5007 lines, presumably a result of 
local jet/cloud interactions in the inner NLR. If jet/cloud 
interactions are responsible for producing the double- 
peaked narrow emission lines, then CXOJ1426+35 is 
likely to b e a strong radio source, similar to the radio 
sources in iHolt et aD ([20Q8.) , and would therefore have 
also been detected. However, the lack of a detection 
makes the presence of (at least strong) jet/cloud interac- 
tions highly unlikely. 

3.3. Rotating Accretion Disk 

A rotating disk with a line emitting surface is capable 
of producing double-peaked emission. The relative posi- 
tions of the peaks and the shape of the lines are depen- 
dent on geometrical factors which can enhance the effect, 
such as the orientation of the disk plane relative to the 
observer's line of sight, and the disk outer radius. Since 
this disk model has been successfully applied to accre- 
tion disks around SMBHs, providing a possible explana- 
tion for the asymmetric (non- Gaussian), broad and often 
double-peaked Balmer emission lines of some AGN and 
quasars, including some candidate binary SMBHs (see 
Section [T]), we consider this scenario for CXOJ1426-I-35 
but find that it is unlikely: the 2-D spectroscopy re- 
veals that the line peaks are at a projected separation 
of ^ 5.5 kpc. This is much larger than the sub-pc sizes 
of such disks based on S MBH masses and Ha line profiles 
(jChen fc Halpernlll989( ) and makes the single accretion 
disk interpretation of CXOJ1426-I-35 unphysical. 

3.4. Unusual NLR Geometry 

One may imagine a disk-shaped NLR which can, in 
principle, produce the double-peaked lines in a way anal- 
ogous to the rotating disk in Section 13.31 One advan- 
tage of this scenario is that NLRs are much larger than 
accretion disks and could explain the observed separa- 
tion. Wh ile a size of 5.5 kpc is la rge for a typical Seyfert 
galaxy (jSchmitt fc Kinnevlll996l ). luminous quasa rs have 
been observed with NLRs of up to 10 kpc (Bcnncr t et al.l 
[200l . The total [OHI] 5007 luminosity of CXOJ1426-f 35 
is nearly an order of magnitude l arger than tha t of th e 
brightest quasar in the sample of iBennert et al.l (|2002| ). 
suggesting that the NLR of CXOJ1426-I-35 may plausi- 
bly extend to the observed separation of 5.5 kpc. 

However, an additional and important consideration 
involves line ratios and line strengths: in this scenario 
both line systems should have similar line ratios and line 
strengths since the ionizing continuum incident on the 
gas and gas densities would be the same since they would 
be part of the same NLR. This scenario was found to be 
plausible for a subset of SDSS double-peaked AGN with 
even [OHI] fiux ratios in ISmith et al.l ()2011| ). However, 
from Table [21 it is clear that the line ratios and line 
strengths, particularly for CIV1549 and [OHI] 5007, are 



not similar. Unless one side of the NLR suffers more ex- 
tinction than the other, and there is no reason to expect 
this, the different line strengths are inconsistent with 
both line systems originating from the same NLR. 

3.5. Accretion Disk Outflows 

Winds moving in opposite directions and ionized by 
a central continuum source can, in theory, produce the 
double-peaked pr ofile of the narrow emission lines (see 
the d iscussions in Xu &: Komo ssa 2009 and Fi scher et al.l 
120111 ). Mechanisms for producing such an out flow in- 
clude radiation pressure acting on gas and dust ([Everettl 
[200l or a hot wind that entrains the NLR clouds 
(jEverett k Murravl l2007| ). A potential problem with 
this scenario is the difficulty of the winds extending 
to typical NLR scales in order to affect the whole 
NLR. However, a SMBH accreting at a relatively high 
Eddington ratio (/Edd = ^BOLZ-^^Edd) may be capa- 
ble of producing a wind strong enough to extend to 
kpc scales. In general, luminous. Type 1 unobscured 
quasars at z > 0.5 may have typical E ddington ratios 
of fadH » 0.25 (iKollmeier et all 120061) ■ though recent 
work bv lKellv et al.l ()2010[ ) suggests that /Edd peaks near 
0.05, or alternatively AGN might instead accrete close 
to Ed dington in a 'buried' phase as described in iKingl 
()2010l). In the sample of nine narrow line Seyfert 1 galax- 
ies (NLSls) with [OHI] offsets in Komossa et al. (2008i) 
the Eddington ratios are high (/Edd = 0.5 — 1.5, with 
the average value /Edd = 0.9), and the authors remark 
that these high Eddington ratios may capable of produc- 
ing the strong outfiows needed to explain the observed 
velocity shifts. In principle, a very rough estimate of the 
Eddington ratio can be made using the bolometric lumi- 
nosity estimated in Section !^. 1. H and estimating a SMBH 
mass assuming all of the stellar mass (calculated in Sec- 
tion l2.3p is within the bulge an d using the Mbh -^Buigo 
relation (jHaring fc Rixl |2004| ) . However, the resulting 
mass of 1.2 x 10^ Mq, and the relation LEdd = 1-3 x 10^* 
Mbh/M0 imply that the bolometric luminosity is ^ 50 
times that of the Eddington limit. This result seems un- 
physical, possibly because this galaxy is still in the pro- 
cess of forming and the Mbh ~ M^uigc relation has yet to 
be established in this system. However, a different test of 
the outflow scenario uses the optical/UV emission lines 
as diagnostic indicators of the ionization parameters and 
kinematics of the 'blue' and 'red' systems. We perform 
two calculations based upon the emission line parameters 
that provide clues as to the origin of the double-peaked 
lines. 

Ionizing Continuum: We can use the ionizing contin- 
uum to estimate the maximum distance at which the 
emission line regions can be so ionized to determine if the 
outflow scenario is plausible. The ionizing continuum is 
described by Q[H'^], which we have estimated from the 
total H/3 luminosity of both the 'blue' and 'red' sources 
via the relation Lhi3 = hi^Hi3{ct^p/aB){^/4:T!')Q[H^]- 
The expected distance can then be found from the re- 
lation U = Q[H'^]/{4:7rr^cnH) where U is the ionization 
parameter estimated as a functi on of gas density from 
the [OII]3727/[OIII]5007 ratios (Komos sa et"aI1 l200l . 
The calculated distances are highly dependent on the 
global covering factor, fl, of the winds as viewed from 
the centr al source. Therefore, w e adopt a lower limit of 
n = 1.6 (|Blustin fc FabianI [20091 ) for each source, which 
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corresponds to an upper limit on the distances. As- 
suming a pure Hydrogen nebula and complete ionization 
{riH = rif,), we estimate the distances for a typical range 
of NLR electron densities {ue — 10^ — 10^ cm~^), where 
10"^ cm^'^ is an upper limit based on the critical den- 
sity for de-excitation of [Oil] 3727, which is well-detected 
in the spectrum of CXOJ1426+35. Given the range in 
electron densities, we find maximum distances of 2 — 5 
kpc and 2 — 6 kpc for which the 'blue' and 'red' systems, 
respectively, can be so ionized. While these distances 
are consistent with the projected separation of 5.5 kpc, 
making the outflow scenario plausible on these grounds, 
it should be noted that these estimates are intended as 
upper limits, and that the projected separation is a lower 
limit on the actual physical separation. Therefore, these 
calculations merely show that within the uncertainties 
of the geometry of the system, two clouds with the ob- 
served [Oil] 3727/ [OIII] 5007 line ratios and H/3 luminosi- 
ties could be ionized by a single source. The limits of 
this argument are apparent when one considers the sig- 
nificantly different line ratios and line strengths between 
the two systems as noted in Section 13.41 In particu- 
lar, the significantly different CIV1549 and [OIII]5007 
strengths are not expected if the outflowing clouds are 
at the same distance from the AGN and have similar 
compositions. However, in principle, a clumpy obscuring 
geometry could account for this observation. 

Ionization Stratification: The outflow scenario should 
produce a positive correlation between velocity offset, 
AY, and ionization potential, LP., otherwise known as 
an 'ionization stratification.' This is because the high 
ionization lines will originate closer to the central source 
where the flux is large enough to produce such lines, 
whereas low-ionization lines will originate further out 
where they are less affected by the outflow. Such a 
line stratification is evident in th e sample of NLSls wit h 
blueshifted [OIII] 5007 lines from KomossaeLahl (|2008aD . 
We investigate AV vs LP. for CXOJI426-h35 using the 
emission lines for which the central wavelength was al- 
lowed to vary during the fitting procedure. From a Spear- 
man rank correlation test, we find probabilities for a cor- 
relation between AV and LP. of 55% and 60% for the 
'blue' and 'red' systems, respectively. These probabili- 
ties indicate that a significant correlation is not appar- 
ent given the available line measurements. In Figure |9] 
we plot AV vs LP. for each system along with the lin- 
ear least squares best fit. In both systems, the velocities 
are all consistent within their errors, with the exception 
of CHI] 1909. Since the CHI] 1909 apparent blueshift ap- 
pears in both systems it may be related to the anoma- 
lous CHI] 1909 blueshift seen in some quasars as a result 
of the increased intensity of the A1907 transistion rel- 
ative to A1909 300 km s^^ offset) at low densities 
(jFerlandl[l981l) . While there is no appar ent cor r elation , 
it should be noted that in the Komoss a et al.l ()2008a[ ) 
sample, two galaxies have the [FeX] line from the coro- 
nal line region (I.P.= 235 eV), significantly higher than 
the highest LP. line observed for CXOJI426-H35 ([NeV], 
I.P.= 97 eV). For the four NLSls in that sample that do 
not include the [FeX] line, the trends are comparable to 
that of CXOJ1426+35. This indicates there is a possibil- 
ity that a line stratification might become apparent with 
a wider range of I.P.s. 



3.6. Galaxy Merger and Dual SMBHs 

A system of two SMBHs within the host galaxy, with 
at least one of them active, would naturally explain the 
two emission line systems as each SMBH would provide 
its own NLR. Since such a system would arise from a 
galaxy merger, one may speculate that it is more likely 
to be found in a disturbed system. Indeed, the morphol- 
ogy of CXOJ1426-I-35, particularly the "lumpy" struc- 
ture evident in the NIRSPEC image in Figure [3 sug- 
gests a possible disturbance of some kind. This possibil- 
ity seems even more likely when one considers that, at 
z = 1.175, galaxy mergers were more frequent than in 
the local universe (|Bridge et al.l[2010l ). We also specu- 
late that the X-ray spectrum (which implies a high level 
of dust extinction) is indicative of a galaxy that has re- 
cently undergone a merger. Therefore, the interpretation 
of CXOJ1426+35 hosting two SMBHs is quite plausible 
and worthy of consideration. We consider two possibil- 
ities involving a galaxy merger and formation of a dual 
SMBH system. 

One AGN Illuminating Two NLRs: This scenario 
was also sugge sted for a candidate dual SMBH in 
iXu fc Komossal (|2Q09f ). It is possible that after the 
merger, only one of the SMBHs would be active but both 
would have their own NLR. The active AGN may there- 
fore be able to ionize both NLRs which would have a 
velocity difference imparted to them by the merger. For 
the NLR around the inactive SMBH, assuming it is com- 
pletely ionized by the active SMBH and with — 10^ 
cm~^, we estimate a free-fall timescale of ~ 5 x 10^ 
years. This is a relatively short time, and it is un- 
likely that we would see the system during this win- 
dow. However, a more rigorous test for this scenario uses 
the expectation that the ionization parameter and Hy- 
drogen number densities should be consistent with one 
NLR being further from the ionizing source than the 
other. In Section l375l we estimated a range of distances 
based on the ionization parameters (estimated from the 
[OII]3727/[OIII]5007 ratio) and a typical range of elec- 
tron densities [ue = 10^ — 10'' cm~^). It is theoreti- 
cally possible to obtain an independent estimate for the 
electron densities from the fiux ratios of emission lines 
produced by different transitions within the same ion, 
particularly [OIII](jA4959 + jA5on7)/jA4363 and assumin g 
a temperature of the gas (jOsterbrock fc Ferlandl l200"6l) . 
However, while the measured [OIII] ratios are well-within 
the typical values for AGN NLRs, we note that is 
strongly dependent on the NLR temperature, and pho- 
toionization models imply NLRs may have multiple zone s 
with different temperatures and densities (|Kroliklll999D . 
Therefore, we have chosen to simply use the typical range 
of densities in our estimates. Given the range of distances 
for each system, the possible ratios between the distance 
of the 'blue' component to that of the 'red' component 
range from 1:3 to 2.5:1. In no case does the ratio of dis- 
tances imply one source is significantly further than the 
other, providing no strong evidence for just one source 
being associated with an AGN. 

Two AGNs: It is generally agreed that galaxy merg- 
ers, through tidal interactions, are effective at enhancing 
SMBH activity by tunneling gas and dust t oward the 
central region of the new d istribution of mass ()Hernquisti 
[1981 ISpringel et al.l[2005l) . Therefore, ff CXOJ1426-h35 
did experience a wet merger, then the merger would have 
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increased the probability for the two SMBHs to be active 
as gas would be forced toward the potential well of each 
SMBH. On these grounds a dual AGN scenario is highly 
plausible and would naturally explain the two emission 
line systems. A potential piece of evidence against this 
hypothesis is the NIRC2-A0 image which reveals only 
one source at rest-frame l/im as opposed to the two nu- 
clei expected in the dual AGN scenario. However, there 
are several properties of CXOJ1426-f 35 which point to a 
high level of internal extinction capable of obscuring an 
AGN continuum at rest-frame 1/xm. In particular, the 
high X-ray HR implies there is at least one AGN suffering 
significant absorption, e.g. nn = 6 x 10^'^ cm^^ (Section 
12.1. ip and modeling of the SED (Figure ^ reveals that 
an AGN with comparable absorption would be virtually 
invisible at rest-frame l/xm (Section 12. 3p . This AGN 
would dominate in the X-ray and IRAC bands, while 
a second, less obscured AGN would be associated with 
the source seen at 1/xm. We note t hat another candidate 
dual AGN (jComerford et al.1l2011[ ) also shows only a sin- 
gle component in AO imaging, but which may host two 
heavily obscured (Compton thick) AGN as revealed from 
Chandra data. As discussed in Section [1] s everal con- 
firmed dual AGN ( Guainazzi ct al. 2005; Bi anchi et al.l 
l2008HPiconcem et al.li2010i:iKoss et al.li2011i) provide ex- 
amples of this scenario, showing only one visible nu- 
cleus, with the other heavily obscured and identifiable 
only through X-ray observations. Given the prevalence 
of heavy obscuration in other candidate and confirmed 
dual AGN, the lack of detection of two cores would not be 
surprising if CXOJ1426+35 in fact hosted a dual AGN. 

3.7. Multiple Scenarios 

Finally, we note that more than one of these phys- 
ical mechanisms may be present. In particular, we 
note the possible coexistence of a galaxy merger and 
AGN outfiows. This follows from the morphology of 
CXOJ1426-I-35, which suggests a possible merger, and 
from the notion that galaxy mergers are potential insti- 
gators of enhanced SMBH activity. If the merger has re- 
sulted in an AGN with a relatively high Eddington ratio 
(e.g., > 0.25 and perhaps close to the Eddington limit), 
the accretion disk radiation pressure may be capable of 
produci ng t he double-peaked narrow lines as discussed in 
Section 13.51 While an ionization stratification is not ap- 
parent given the current data, we could not rule one out 
for either system. This scenario is consistent with both 
the morphological suggestion of a recent merger while al- 
lowing for the existence of strong outflows driven by the 
radiation pressure of AGN accretion disks. 

3.8. Summary 

We have discussed six separate physical mechanisms 
that are each capable of producing double-peaked emis- 
sion lines, and we have discussed the likelihood of each 
scenario based on the available data. We find that a 
superposition, even within a galaxy cluster, is highly un- 
likely given the small physical and velocity separation, 
and that the rotating accretion disk scenario is highly 
unlikely given the several kiloparsec separation. We find 
that a jet/cloud interaction is also highly unlikely sim- 
ply because it is not a strong radio source. The scenario 
of a NLR with an disk-shaped component can naturally 
explain the double-peaks, but the disagreement between 



the line ratios and line strengths of both systems is not 
expected in that scenario. We find that the biconical 
outfiow from an accretion disk can not be ruled out, and 
that the primary limitation is the power of the accretion 
disk. A recent galactic merger would potentially make 
the biconical outfiow a more likely explanation, although 
there is no clear line stratification. In the case of a single 
AGN, the SED modeling implies there might also be a 
very young stellar population present which accounts for 
the UV excess. A merger would also increase the possibil- 
ity of a single AGN illuminating the NLR of each galaxy, 
though this scenario may not be likely given the short 
free-fall timescale and the similar distances of each cloud 
based on a single ionizing source. A more likely scenario 
is that there are actually two AGN present. The lack of 
two nuclei in the NIR AO image provides possible evi- 
dence against a second AGN, but this may be naturally 
explained by the high obscuration, inferred from both the 
X-ray data and the SED modeling, which would obscure 
the continuum radiation. This scenario is also consistent 
with the several kiloparsec separation of the AGN com- 
ponents, since many of the dual AGNs already identified 
have comparable separations, and it is consistent with 
the AGN-like line ratios. 

4. TESTS FOR DUAL AND BINARY SMBHS 

Based upon limited data, determining the true nature 
of many candidate dual and binary SMBH systems has 
proven to be difficult given the attractiveness of alterna- 
tive explanations. Among the current handful of candi- 
dates there is a variety of data and interpretations, many 
of which have been tested and often with ambiguous re- 
sults. Here we use CXOJ1426-I-35, in comparison with 
other candidates, to review how these candidates can be 
tested. 

4.1. One- Dimensional Spectroscopy 

Given the availability of spectroscopic databases that 
include large samples of AGN and quasars, emission line 
diagnostics are commonly used to select the most promis- 
ing double-peaked emission line candidates for follow-up 
observations. From one-dimensional spectra alone (i.e. 
SDSS fiber spectra), much can be done to address the 
likelihood of a dual/binary AGN. Emission line ratios can 
provide informatio n about the source of the ionizing con- 
tinuu m radiation (jBaldwin et al.|[T981l : iKaufFmann et al] 
l2003f ) , and careful examination of spectra may also reveal 
asymmetries, velocity-offsets, or multiple components, 
thereby providing valuable information about the kine- 
matics of the NLR gas. We have performed these analy- 
ses for CXOJ1426-I-35 in the above discussion, and they 
have been done for all of the other candidate dual/binary 
SMBHs yet identified, providing valuable information 
and resulting in a variety of interpretations. For example, 
in J0927-I-2943, close examination of the optical spectra 
(which features two prominent emission line systems, one 
with just narrow lines and the other with broad and nar- 
row lines) has allowed for multiple mo dels to be put forth, 
including a r ecoiling SMBH (iKomossa et al.ll20081: ) , a b i- 
nary SMBH (Bogdanovic et al. 2009b: iDotti et al.ll2009D . 
a superposition (Shields et aL 20Q9ai ) and a large and 
smal l galaxy interac t ing n ear the center of a rich clus- 
ter (iHeckman et al.l I2009D. As another example, in 
J1316-H753 lXu fc Komossal (|2009h identified two NLRs 
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and at an intermediate redshift there is a broad com- 
ponent to each emission hne, which the authors discuss 
may be evidence of a starburst-driven outflow that was 
triggered by a merger. 

We also note that spectroscopic observations over mul- 
tiple epochs could potentially reveal any orbital motion 
in a close binary SMBH system. Though CXOJ1426-I-35 
does not fall into this class as the implied separation 
is much too large, candidate close binary SMBHs could 
be identified in this manner. This test ha s been per- 
formed for b oth J0927-H2943 (iVivek et all [20091 ) and 
J1536-f0441 (jChornock et all |2010D which were iden- 
tified as candidate sub-parsec binari es with separa- 
tions of ~ 0.34 pc (iDotti et all [20091) and - 0.1 pc 
([Boroson fc Lauen 12009). respectively. In both cases, 
the acceleration is found to be zero within the errors. 
However, it has been pointed out that such candidates 
would be preferentially identified when their rad ial ac- 
celerations are smallest ([Bogdanovic et al.1l2009al ). 

4.2. Two-Dimensional Spectroscopy 

By providing information about the spatial separation 
of the two emission line components, two-dimensional 
spectroscopy has proven to be highly effective at identi- 
fying the most promising dual AGN candidates from the 
much larger sample of double-peaked narrow emission 
line sources. In this way, CXOJ1426-I-35 has been shown 
to be a strong dual AGN candidate based on the two 
emission line components spatially separated and with a 
velocity difference of several hundred km s~^. At lower 
redshifts, there are many other similar candidates for 
which two-dimensional spect roscopy reveals two clea rly 
separated components (e.g. iComerford et al.l l2009al fbl). 
Two-dimensional spectroscopy ruled out the close binary 
hypothesis for J0927-I-2943, but show ed that it may sti ll 
be a system of two merging galaxies (jVivek et al.ll2009f ). 
A large sample of double-peaked AGN with follow-up 
two-dimensional spectroscopy will increase the number 
of candidate dual AGN (Comerford et al., in prep). How- 
ever, as made evident in our analysis of CXOJ1426+35, 
two-dimensional spectroscopy alone can not confirm a 
dual AGN. 

4.3. Morphology 

Morphology can provide powerful evidence for or 
against the dual AGN hypothesis. A number of ob- 
jects which show double-peaked emission lines have 
been imaged in optical and near-IR bands, where the 
results have shown morphological evidence of inter- 
actio ns /mergers such as tida l tails and multiple nu- 
clei (lJunkkarine n et all 120011; I Comerford et al.l l2009bl : 
iGreen et al.l 12010: Liu et al. 2010a). From HST imag- 
ing, the galaxy COSMOS J100043.15-H020637.2 h as been 
establ ished as hosting eithe r a dual dC omcrfor d et al.l 
[2009a) or a recoihng SMBH (Civano et al. 2010). In th e 
case of EGSD2 J142033.66-i-525917.5 (Gerke et al.. i20?)7l ). 
an HST ACS image shows morphological evidence of a 
possible inter action, but no cle ar evidence f or mu ltiple 
nuclei. Rosari o et all (j20lH ) and lShen et al.l ()201lD have 
obtained NIR imaging of several [OIII] double-peaked 
AGN where the results have shown a variety of possible 
origins of the double-peaked narrow emission lines, in- 
cluding some sources with morphological evidence sug- 
gesting the most likely interpretation is a dual AGN. 



iFu et al.l ()2010[ ) present images for a sample of several 
such candidates which reveal multiple nuclei in some 
cases, though they caution that spatially resolved spec- 
troscopy is_j;e2uired_to confir m any case as a true dual 
AGN. In lMcGurk et al.l (|2011l) integral field spectroscopy 
was able to confirm a dua l AGN. That so urce was con- 
firmed in the same way in iFu et al.l ()2011[ ). though their 
sample revealed that, out of 42 sources, in all but two 
the double-peaked lines were produced by gas kinemat- 
ics. These studies have highlighted the usefulness of 
integral field spectroscopy in eliminating the ambigui- 
ties left by double-peaked emission lines and even high- 
resolution imaging. In the case of CXOJ1426+35, NIR 
AO imaging revealed only a single nucleus at a rest-frame 
effective wavelength of 1/im, providing possible evidence 
against the dual AGN scenario. As with other dual AGN 
candidates, integral field spectroscopy will be useful for 
CXOJ1426-F35. 

In the case of heavy nuclear obscuration, however, de- 
tection of double point sources in X-ray images of close 
ULIRGS has bee n effective at identif ying dual AGN, 
e.g., NGC 6240 (IKomossa et al.l [200l. ESO509 -IG066 
(iGuainazzi enill I2005D 3C 75 iH udson et al.l [2006) . 
Mrk 436 (iBianchi et al.1 [20081). and IRAS 20210-M121 
(jPiconcelli et al.l 20101) . The advantage of this method of 



detection is that hard X-rays are much less attenuated 
by intervening material in the galaxy, and recent mergers 
are known to have high column densities in the nuclear 
regions. This might be the case for CXOJ1426-I-35, as 
dust may have a strong role in hiding an AGN near 1/xm, 
whereas two AGN would be revealed from hard X-rays. 
Though CXOJ1426+35 is unresolved in our image, this 
is due to the off-axis detection and the resolution of the 
ACIS detector. With the HRC, the two X-ray compo- 
nents could potentially be resolved in CXOJ1426-I-35. 

5. CONCLUSIONS 

We have analyzed all available data, new and archival, 
of the source CXOJ1426-t-35 which has an intriguing 
spectrum showing double-peaked optical/UV emission 
lines. The two emission line components, separated spa- 
tially by 0.69" (5.5 kpc) and in velocity-space by 690 
km s~^, make CXOJ1426-I-35 similar to other candidate 
dual AGN, though it is unique in that it is the highest- 
redshift candidate dual AGN yet identified. The opti- 
cal/NIR imaging shows an elongated and "lumpy" mor- 
phological structure suggestive of a disturbance, whereas 
the AO image shows only a single source. We have 
discussed several interpretations of the physical nature 
of CXOJ1426-I-35, including a chance superposition, a 
jet/cloud interaction, line emission from a rotating ac- 
cretion disk, an unusual NLR geometry, a strong bicon- 
ical outflow from an accretion disk, a dual SMBH and 
a combination of these scenarios. We find that the only 
plausible scenarios are a strong outflow, a dual SMBH 
system, or a combination thereof. If CXOJ1426-I-35 re- 
cently underwent a merger, an increased Eddington ratio 
may have resulted, making the disk outflow scenario a 
possible explanation. Alternatively, a merger may have 
resulted in the actual case of a dual AGN, or, less likely, 
a single AGN iUuminating the NLRs of both SMBHs. 
Regardless, CXOJ1426-)-35 is most likely a merger rem- 
nant. 

Finally, we compared CXOJ1426-I-35 to other candi- 
dates with the aim of reviewing methods of testing the 
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dual AGN hypothesis. Analysis of single-epoch ID and 
2D spectra provides crucial information about the ion- 
izing continuum and the velocity and spatial separation 
of the two emission line components, but high-resolution 
imaging and multi-epoch spectra have proven most ef- 
fective at confirming or rejecting the dual/binary SMBH 
hypothesis. GXOJ1426-I-35 represents an interesting ex- 
ample of a dual AGN candidate, with many proper- 
ties similar to other candidates in the literature, but 
also with several unique characteristics. In particular, 
CXOJ1426-K35 is the highest redshift dual AGN candi- 
date to-date, at a redshift of z = 1.175 when galaxy 
mergers were more frequent, making the dual AGN hy- 
pothesis especially appealing. This makes it a fascinat- 
ing addition and will help in developing a better under- 
standing of how to look for the potential signatures of 
dual/binary SMBHs and of the often complex emission 
line regions of AGN. 
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Fig. 1. — Chandra/ ACIS spectrum and model residuals of CXOJ1426+35 in the observed frame. The top panel shows the combined data 
sets of the extractions from both the 120 ks and 58 ks images (in Wang ct al. 2004 the 58 ks image was corrected to 52 ks, hence their 
quoted total time of 172 ks), and the best fit absorbed power-law model is overlaid. The bottom panel shows the residuals from the best 
fit model. 




Fig. 2. — Stacked g -\- r i SDSS image centered on %ic position of 'fcxOJ1426+3^.° The larger image is 90" on a side with North up and 
East to the left. The inset is a close-up of CXOJ1426+35, 10" on a side, showing that the galaxy is clearly extended with a position angle 
of P.A.= —54.1°. Boxes show the spectroscopic slit configurations obtained using LRIS (P.A.= —21°, slit width = 1.5") and NIRSPEC 
(P.A.= -54.1°, slit width = 0.7"). Note the bright star, r = 13.04, ~ 25" to the east of CXOJ1426-I-35, making this galaxy an ideal AO 
target. 
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Fig. 3. — Prominent emission lines seen in the Keck/LRIS 2-D spectrum of CXOJ1426-I-35. The two components arc most clearly 
seen in the forbidden lines which arc narrower and less blended. For lines with blended doublets, particularly CiV1549 for which 
the 'blue' component is much weaker than the 'red', the two components are more difficult to resolve but still present. From left 
to right, the first four panels show emission lines from the blue side of the dichroic, and the last five panels are from the red side 
of the dichroic. The dispersion scale is 1.90 A pixcl^-*^ for the blue side and 1.18 A pixel^^ for the red side. Note that the blue 
side was binned by two pixels in the dispersion direction. Along the spatial axis the scale is 0.135 arcsec pixel^^. In each box the 
horizontal axis (relative velocity) spans 4000 km s"'^ and the vertical axis (spatial separation along the slit) spans 10". At a P. A. 
of —21°, the spatial separation between the two components is 0.48" and the velocity separation is 790 km s~^, as measured from [NeV]3426. 
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Fig. 4. — Keck/LRIS spectrum of CXOJ1426+35, in the observed frame, extracted from a 3" diameter region including both velocity 
components, mimicking the source spectrum through an SDSS fiber. The emission lines and atmospheric absorption lines (A and B band) 
are labeled. Note the double high ionization narrow emission lines, suggestive of two Type 2 AGN in the process of merging. Two 
dimensional spectra of these emission lines are presented in Figure [S] 
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Fig. 5.— From left to right: stacked g + r + i SDSS (0.40 arcsec pixel"!), Y-band NIRSPEC (0.18 arcsec pixel"!) ^^^^ i^'-band NIRC2 
- AO (0.01 arcsec pixel"!) jj-Q^ges. In each image, the scale bar represents 10 kpc at ^ = 1.175. The AO image has been smoothed by 
convolving with a Gaussian function using the 'gauss' task in IRAF. 
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Fig. 6. — Panel showing the H/3 and [OIII]5007 emission lines in the 2-D NIRSPEC spectrum taken at P.A. of -54.1°. The dispersion 
scale is 2.10 A pixel"^, and the scale along the spatial axis is 0.143 arcsec pixel"^. In each box the horizontal axis (relative velocity) 
spans 3000 km s~^ and the vertical axis (spatial separation along the slit) spans 9". At P.A. of -54.1° (aligned parallel to the major axis 
of the galaxy), the spatial separation between the two components is 0.69" and the velocity separation is 690 km s~^, as measured from 
[OIII]5007. 
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Fig. 7. — NIRSPEC spectrum of CXOJ1426-I-35, in the observed frame, extracted from an aperture 3" in diameter that includes both 
velocity components. The emission lines are H/3, [OIII]4959 and [OIII]5007. Overlaid on the spectrum are the best-fit Gaussians. Each 
emission line is fit by two Gaussians, one for the 'blue' and one for the 'red' emission line component. 
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Fig. 8. — SED fits to the broad-band photometry assuming a single {panel a) and dual {panel b) AGN components. The best-fit AGN 
components in each case are shown by the dotj-dashed lines, while the Im component is shown by the dashed lines and the Sbc component is 
shown by the dotted lines (sec Asscf ct al. 2010^, for details on the SED templates). The solid line shows the best-fit SED, which corresponds 
to the addition of the different components. 
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Fig. 9. — Plot of emission line velocity-splitt ings fo r the 'blue' and 'red' systems relative to 2 = 1.1751 (i.e. the systemic velocity as 
measured from the H9 absorption line in Scction [2.2.1|l . where positive velocities indicate a blueshift. The plotted data points correspond to 
emission lines for which the positions were allowed to vary during the fit. Some of the line labels have been shifted for clarity. The straight 
lines are the linear least squares best fits. The quoted probabilities for each system are the two-sided significances, based a Spearman rank 
correlation, of the deviation from the null hypothesis that there is no evolution of AV with LP. In both systems, the probabilities provide 
no indication of a corr elation. The apparent blueshift of GUI] 1909 in each system may be the result of an increase in the 7(A1907)//(A1909) 
ratio (see Section 13.511 . 
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TABLE 1 

Photometric data for CXOJ1426+35. 
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22.49l«°^ 


SDSS 
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5975 A 
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26.16^21? 


NOAO/NEWFIRM 
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25.34t«-;}? 
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oo.oo_j^4 23 


NOAO/NEWFIRM 




1.7^m 
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45.77ti;l 


Spifaer/IRAC 
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16, 


.21±0, 
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Spitzer/IRAC 
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Magnitude determined for the 3.8" radius aperture corrected to total magnitude (see Section l2. 1.211 . 
b SDSS Model magnitudes (see Section [2X3t . 

Magnitudes are for a Kron-like elliptical aperture (see Section l2. 1.31 1. 
^ Magnitudes determined for the 4" radius aperture corrected to total magnitudes (see Section l2.1.4|l . 
° The Chandra flux densities are not absorption-corrected. 



TABLE 2 

Fluxes, FWHMs, redshifts and velocity-splittings for emission lines in the both systems. 



Line 


Flux 


Blue System 
FWHM 


2 


Flux 


Red System 
FWHM 


z 


Velocity-splits 


CIV1549 


23.30±6.23 


1120.95'' 


1.1718^^ 


109.04±9.85 


1120.95±85.65 


1.1776 


797.37±75.49 


HeII1640 


20.70±6.50 


818.90±162.03 


1.1715 


29.52±6.58 


872.33±141.95 


1.1779 


882.47±157.33 


[OIII]1663 


1.65±0.47 


555. 68^ 


1.1720^ 


3.10±0.54 


474. 17^ 


1.1770'' 


690.95'' 


cm] 1909 


14.72±6.19 


925.58±154.08 


1.1687 


14.34±6.22 


1003.31±192.97 


1.1745 


971.68±236.55 


CII]2326 


3.30±1.90 


880.83±63.68 


1.1718 


2.32±1.98 


901.53±239.35 


1.7758 


799.66±462.92 


[NeV]3426 


19.14±1.25 


590.33±39.64 


1.1715 


17.48±1.20 


509.30±36 


1.1773 


795.14±28.90 


[OII]3727 


178.05±76.04 


719.33±170.34 


1.1722 


80.10±82.79 


697.81±237.48 


1.1771 


681.88±333.59 


[NeIII]3869 


67.97±5.73 


615.60±63.09 


1.1718 


37.06±6.54 


542.78±97.99 


1.1772 


736.66±49.73 


[NeIII]3968 


36.13±2.27 


615. 60^^ 


1.1718'= 


19.35±1.61 


542.78'= 


1.1772'= 


736.66" 


H7 


59.11±8.21 


587.03<i 


1.1724^ 


33.46±6.39 


474.17'' 


1.1766'' 


579.46'' 


[OIII]4363 


26.47±5.04 


555.68'' 


1.1720'' 


12.34±3.46 


474.17'' 


1.1770'' 


691.04^ 


H/3 


178.67±63.15 


587.03±274.57 


1.1724 


< 79.98 


474.17'' 


1.1766 


579.46±156.43 


[OIII]4959 


511.11±108.30 


487.90±59.76 


1.1721 


186.61±40.34 


474.17'' 


1.1770'' 


675.94±44.75 


[OIII]5007 


1581.48±132.82 


: 555.68±19 


1.1720 


477.93 


474.17±45.55 


1.1770 


690.95±44.82 



Note. — Fluxes arc in units of 10 '^ erg s ' cm ^ A ', and FWHMs and velocity-spittings relative to z = 1.175 are 
in units of km s~'. All errors correspond to Icr uncertainties. 
^ Fixed to the 'red' CIV1549 component. 
'' Fixed to [OIII]5007. 
" Fixed to the value for [NeIII]3869. 
'' Fixed to the value for H/3. 



